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As a unifying framework for examining several properties that nominally involve eigenval¬ 
ues, we present a particular structure of the singular values of the Gaussian orthogonal 
ensemble (GOE): the even-location singular values are distributed as the positive eigen¬ 
values of a Gaussian ensemble with chiral unitary symmetry, while the odd-location 
singular values, conditioned on the even-location ones, can be algebraically transformed 
into a set of independent %-distributed random variables. We discuss three applications 
of this structure: first, there is a pair of bidiagonal square matrices, whose singular 
values are jointly distributed as the even- and odd-location ones of the GOE; second, 
the magnitude of the determinant of the GOE is distributed as a product of simple 
independent random variables; third, on symmetric intervals, the gap probabilities of the 
GOE can be expressed in terms of the Laguerre unitary ensemble. We work specifically 
with matrices of finite order, but by passing to a large matrix limit, we also obtain new 
insight into asymptotic properties such as the central limit theorem of the determinant or 
the gap probabilities in the bulk-scaling limit. The analysis in this paper avoids much of 
the technical machinery (e.g. Pfaffians, skew-orthogonal polynomials, martingales, Meijer 
G-function, etc.) that was previously used to analyze some of the applications. 
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1. Introduction 

This paper studies the structure of the singular values of the Gaussian orthogonal 
ensemble (GOE), using it as a unifying framework for examining several properties 
that nominally involve eigenvalues. Here, the GOE n of order n is the ensemble of 
real symmetric random matrices 

G = (X + X')/2, 

where X is an n x n Gaussian matrix with all entries independent standard normals. 
Since the singular values of symmetric matrices are the magnitudes of the eigenvalues, 
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the ensemble of singular values will be briefly denoted by |GOE„|. Central to our 
discussion is the immediate set decomposition, 

|GOE n | = even |GOE„| U odd|GOE n |, (1.1) 

of the ordered singular values according to the parity of their indices, where the 
even-location decimated ensemble even|GOE„| is defined by taking the 2nd largest, 
4th largest, etc. singular value, and similarly for odd |GOE n |. 

Our first set of main results relates the decomposition (1.1) to the eigenvalues 
of a Gaussian ensemble with chiral, or anti-symmetric, unitary symmetry. Namely, 
with X as above, the ensemble of real skew-symmetric random matrices 

A = {X-X ')/2 

will be called the anti-GUE with its (almost surely) different and positive singular 
values written briefly as aGUE n (if n is odd, there is a surplus singular value zero, 
which is omitted). a Then, the following structure holds. 

Theorem 1.1. Denoting equality of the joint distribution by =, there holds 

even |GOE„| = aGUE„. (1.2) 

We will give two proofs that differ in their handling of the odd-location singular values: 
one (Section 4) by algebraically transforming them to a set of independent positive 
variables, each distributed as \2 and, if n is odd, a surplus Xii the other (Section 7) by 
integrating them out. Both proofs are based on an algebraic factorization (Section 2) 
of the joint density of |GOE„|, where one factor depends only on the even-location 
singular values, the other on the odd-location ones. If we recall the superposition 
representation, see [8, Eq. (2.6)] or [7, Thm. 1], 

|GUE„| = aGUE n U aGUE„ +1 , 

of the singular values of the Gaussian unitary ensemble (GUE), with both ensembles 
on the right drawn independently, Theorem 1.1 immediately implies the following 
remarkable relation between the singular values of GUE and GOE: 

Corollary 1.1. With the ensembles on the right drawn independently, there holds 

|GUE„| =even|GOE„| U even|GOE n+1 |. (1.3) 

The superposition (1.3) bears a striking similarity with a corresponding result for 
the eigenvalue distributions, see [10, Eq. (5.9)] and [11, Eq. (6.14)], namely 

GUE n = even (GOE„ U GOE„ +1 ). 

Our second set of main results (Sections 4/5) sharpens Theorem 1.1 by realizing 
|GOE n | as the singular values of other matrix models that reveal a rich additional 


a In this paper, the Gaussian weights are e P x2 / 2 with /3 = 1 for orthogonal and (3 = 2 for unitary 
symmetry. 
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structure. A first model (Corollary 4.1), initially identified by comparing moments 
of the product of the even singular values to known moments of the determinant 
[1, Eqs. (23) and (24)], is constructed by bordering the skew-symmetric matrix A 
defining the anti-GUE with an independent standard normal vector b £ that is 
to say, the singular values of G and those of 


H = {b A) 


are both distributed as |GOE„|. To the same end, the bordering vector could also 
be chosen as b = r n e i, where r n is a x n -distributed variable, independent of A, and 
e\ denotes the first unit vector. The precise effect of such borderings on the singular 
values of a matrix is studied in the preparatory Section 3. 

Using a technique (Lemma 5.1) that was, in essence, introduced by Dumitriu and 
Forrester [6], this bordered matrix model is finally (Theorems 5.1/5.2) transformed 
into a pair (R even , R odd ) of bidiagonal square matrices, whose singular values are 
jointly distributed as even |GOE n | and odd |GOE n |. Both matrices depend in a very 
simple fashion on a set of independent ^.-distributed random variables. Specifically, 
for n = 2m even (the structure of the odd order case is similar), we get 


and 


/£°dd 


l 

71 


/TCf + ^2m-2 \ 

(2m-1 6m-4 


V 


6 6 


j^even 


1 

71 


/Cl ?2m-2 \ 

£,2m— 1 (,2m —4 

(5 (2 

V (3/ 


where ( 1,(21 ■ ■ ■ 1 ( 2 m are independent random variables, with distributed as \k- 
The singular values of R odd correspond to odd|GOE 2m | and those of R eveD to 
even | GOE 2to |, both drawn from the same ensemble. 

As a striking application (Corollary 5.1) of this new matrix pair model, we 
establish that |detGOE„| can be expressed explicitly as a product of independent 
random variables. Specifically, for n = 2m even, the determinant of M = \/2G 
factors as 

Idet M\ = (i^(l + 2(l m ■ (l ■ (l ■ ■ ■ CL_r, (1.4) 

with independent variables Cfc distributed as \k- A similar factorization holds in the 
odd order case. The form of these variables explains the absence of large prime factors 
in the moments of the determinant, and leads to a new, simple proof (Section 6) of 
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the known central limit theorem for log |detGOE„|, cf. Delaunay and Le Caer [4, 
Section III], Tao and Vu [17, Thm. 4], While the representation (1.4) of |det M\ as 
a product of independent random variables can be found implicitly in the work of 
Delaunay and Le Caer, namely in form of a factorization [4, Eq. (41)] of the Meijer 
G- function representation of the Mellin transform of det M into hypergeometric 
terms, see the discussion of (5.6), Tao and Vu, who approximated the log-determinant 
by a sum of weakly dependent terms, speculated that such a representation would 
not be possible [17, p. 78]. 

Our third set of main results (Section 8) studies the implication of Theorem 1.1 
on the inter-relation of gap probabilities, that is, the probabilities E(k; J) that the 
interval J contains exactly k eigenvalues drawn from a random matrix ensemble. 
Specifically, for order n, we get 

^goe(2^ + h ~ 1; (— s, s)) + Eq OE ( 2k + /z; (—s, s)) = -E"gue(^! (0, s)), 

where /z = 0,1 denotes the parity of n. This formula was previously known only in 
the case /j = 0, see Forrester [8, Eq. (1.14)]. We initially used a heuristic argument, 
see (8.3), to extrapolate the formula to the case /z = 1. A substantial portion of 
the present discussion was derived from attempts to justify the heuristic, after this 
prediction held up under numerical scrutiny. Taking the bulk scaling limit of both 
cases provides a new, simpler proof of a remarkable formula previously obtained by 
Mehta relating the gap probabilities of the GOE and those of the Laguerre unitary 
ensemble (LUE), see (8.4). 

Notation. In contrast to the previous analyses mentioned, where either ensembles 
of odd (e.g., if Pfafhans were used) or of even order (e.g., if Mellin transforms were 
used) have typically presented considerable technical complications, our treatment 
of ensembles of even and odd order is nearly identical. The formulae themselves, 
however, will often depend on the parity /z of the underlying order n and we will, 
throughout this paper, write 

n = 2m +/i (n = 0,1), rh = to + /i, (1.5a) 

that is, 

m = \n/2\ , m = \n/ 2] , fi = |"n/2] — [_77./2j . (1.5b) 

Terms that only appear for n odd will be written with a factor /z in a sum and 
with an exponent n in a product; etc. This way, without suggesting any natural 
interpolation between the cases /i = 0 and p = 1 (with the notable exception of the 
usage of the heuristic duality principle (8.3) that started our work), we simply avoid 
writing out awkward case distinctions. 

2. Joint Density of the Singular Values 

In this section we establish, in two different ways, the joint probability distribution 
of the singular values cr,- = | Ay | of the GOE induced by the corresponding density 
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for eigenvalues as given by the symmetric function 

n 

p(A i, A 2 , ..., A„) = c n e • |A(Ai, A 2 , .. A n )| (2-1) 

fe =1 

with some normalization constant c n and the Vandermonde determinant 

1 1 1 X 

Cl C 2 ' ' ' Cn t— r 

: ; : 

k>j 

cn— 1 c n — 1 Cn— 1 

SI S2 ‘ ‘ ‘ S n / 

We will frequently use that A(£i,... ,£ n ) > 0 if the arguments are increasingly 
ordered, £1 < ■ ■ ■ < £„. 

By symmetry, we can establish the joint density of the singular values by re¬ 
stricting ourselves to the cone of increasingly ordered singular values 

0 < <7i < • • • < a n , (2.2) 

this way parametrizing |GOE„|. To simplify notation and to avoid case distinctions 
between odd and even order n in later parts of the paper, we introduce two further 
sets of coordinates for this cone. Writing, as detailed in (1.5), n = 2 m + y and 
rh = m + y with /i = 0,1, the coordinates 

Xj = <J 2 j -i (j = 1, • • •, rh), yj = cr 2j (j = 1,..., m) (2.3a) 

satisfy the interlacing property 

0 ^ X! < 1/1 < x 2 < y 2 < • • • T"rh ^ Vm 5 (2.3b) 

formally adding the value y m + 1 = 00 if y = 1. With aV and yl denoting the x and y 

vectors with their components taken in the reverse order, so = (#*, Xm-i, ■ ■ ■, Xi) 
and yl = (y m , ym- 1 , ■ ■ •, 2 / 1 ), we define, depending on the parity of n, the coordinates 

(t,s) = (y 4 -,® 4 -) (/z = 0), (<,s) = (x l ,y l ) (y = l), (2.4a) 

satisfying the interlacing property 

t\ ^ Sl ^ t 2 > S2 ^ ^ tm ^ Sm > 0, (2.4b) 

again formally adding the value s TO +i = 0 if y = 1. A large part of the apparent 

dependence on parity is the fact that some results, like Theorem 2.1, have stable 
expressions in terms of the (x,y) coordinates, while others, like Theorem 4.1, are 
stable in the (t, s ) coordinates. Since the mapping from a = (a i,..., a n ) to either 
the pair of coordinates (a:, y) or (t, s) is orthogonal, transforming the density between 
the three sets of coordinates is simply done by inserting new variable names for old 
ones. Note that the s variables parametrize the even-location decimated ensemble 
even |GOE„| while the t -variables do the same for odd |GOE„|. We call them the 
even and odd singular values. 
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Supported on the cone defined by (2.2), the joint probability density of the 
singular values is 

n 

q(a 1: ... ,a n ) = n\ ^ p(eiCTi,---,e n ^n) = c n n\Y\_e~ crk/2 D(a 1 ,... 1 a n ) 
ee{±i}" fc=l 

with 

D(ai,... ,a n ) = ^2 |A(ei«Ji,..,,e n <j n )|. 

ee{±i}" 

To determine the signs of the Vandermonde terms it suffices to discuss the case 
<t i < • * * < a n : we then get, because sign(efct7fc — ejOj) = e*, if k > j, 

n 

sign A(ei<7i,... ,e n o n ) = e k = ejT 1 = e k . 

k>j k= 2 k even 

Hence, by continuity, there holds on all of (2.2) 

D(a 1 ,...,a n ) = ^2 do(e)A(e 1 <r 1 ,...,e n a n ), 9 0 (e) = e k . (2.5) 

eG{±l} ri k even 

The form of do(e) suggests we proceed in terms of the (x,y) coordinates introduced 
in (2.3). With respect to these coordinates, we obtain the following theorem. 

Theorem 2.1. The joint probability density o/|GOE„|, supported on the cone (2.2) 
and expressed in the coordinates (2.3), is given by 

( rh \ / m \ 

ne-^ /2 -A(x?,...,4)J • \ Y[yk e ~ v2k/2 ■ Mvl, ■ ■ ■ ,yl l )J , ( 2 - 6 ) 

where c n is the normalization constant of the GOE -density (2.1). 

Remark 2.1. Despite the fact that the joint density factors on its domain of 
support, it does not reveal an independence between the underlying variables x and 
y. Their dependence is entirely by the interlacing (2.3b). 

We will give two different proofs of the theorem. The first uses the determinantal 
structure of the Vandermonde terms to establish the factorization, while the second 
uses their polynomial structure and their symmetries. It is our consideration that 
the first proof is more straightforward, while the second provides additional insight 
into the structure of the factorization. 


Proof by Determinantal Structure 

We write D(x;y) for (2.5) when expressed in terms of the (x,y) variables (2.3); it is 
convenient to split the sign changes e into e x and e y accordingly and to use 

0 „(e) = 0(6*0, OU") 
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Since determinants are invariant with respect to a simultaneous row and column 
permutation so that the odd columns and rows occur before the even ones, we 
express the Vandermonde terms as 


A(e< 7 i,..., ecr n ) = det 


rf( 


rtf 


by writing the determinant column-wise with 


t f B \ 


Tr^\x) = 


/x+2 


y x V+2n-2j 


Now, we calculate 
D(x;y) = 

9(e v ) det 

c*e{±ir 

e y e{±i} m 

= ^ 9(e v ) det 

e»S{±l} m 


( x i) ••• 


r (A) 


(^m) 


Jrn), 


A m )l 


rtf 


' no m \yi) ■ 

•• A\ym) 

0 e v A m) {y i) ■ 

•• ^i m \ym) 

(A* = 0 , 1 ). 


no m) (yi) ■■ 

■ A\ym) \ 

e v A m) (yi).. 

■ e m^ m) ( 2 /m)/ 

A\yi) ■■ 

• A\ym) \ 

e i 7r i m) (yi) •• 

■ tArHym)) 


2 m det ... 4 ” l) (^)) ' ^2 ®( eV ) 2 det (^ 7 rJ m) (j/i) ... Tr[ m \y m ) 


£ »e{±i}™ 

= 2 m det (^jTQ n \xi) ... 7 r£ m) (a; A )) • 2 m det ^ 7 T^ r " ) (?/i) ... n[ m \y m 2j . 

By noting m + m = n and by expressing the result in terms of Vandermonde 
determinants, we finally get 

D{x; y) = 2 " • A (a;?,..., ■ y 1 ■ ■ ■ y m A (yf, ...,y 2 m ). (2.7) 

This factorization establishes Theorem 2.1. 


Proof by Polynomiality 

We give a second proof of the factorization (2.7) based on the observation that the 
sum in (2.5) defines a polynomial of homogeneous degree at most Q). We identify 
the factors by symmetrizing known vanishings of this polynomial. Extending the 
definition of D polynomially to all real values of its arguments, we get in particular 


D(e iOi,..., e„cr„) = 9 0 {e)D(a 1 ,...,o n ) (e € {±1}”) 


( 2 . 8 ) 



8 F. Bornemann and M. La Croix 


and, inherited from the Vandermonde terms, D is antisymmetric with respect to 
permutations of either the even or odd indices of er since both sets of permutations 
leave the factor 9o(e) invariant. 

Now, if c Tj = tJj+2 for a belonging to the cone (2.2), we have crj = <Jj+i = Cj+2 
and, hence, by the pigeonhole principle, for each choice of signs e at least one of 

e j a j = e i+i°j+i or e j a j ~ e j+2<Tj+2 or e j+1 a j+ i = Cj +2 a j+ 2 

holds. Therefore, each of the Vandermonde terms in (2.5) vanishes. It follows that 
<jj — (Jj+2 and by (2.8) also <jj + <Jj+2 divide D for all j, thus so does the product 
<jj — a | +2 . We also note that if cr 2 = 0, we have o\ = <r 2 = 0 and, hence, for each 
choice of signs eicri = e 2 cr 2 . Once more each of the Vandermonde terms in (2.5) 
vanishes and it follows that er 2 divides D. 

In terms of the {x, ^-coordinates (2.3a) we thus see that x 2 — x 2 +1 , y? — j/? +1 
and yi divide D(x;y). Invoking the antisymmetry with respect to either x or y, we 
see that D is divisible by — xf., yj — y\ for every j ^ k and by y 3 for every j. 
These factors contribute homogeneous degree 



so D cannot have any other non-unit factors and we get 

D(x; y) = d n ■ A(z?,..., x%) ■ yi • • • y m A(yf,..., y 2 m ) (2.9) 

with a positive constant d n . This is (2.7) except for identifying d n =2™. 

Remark 2.2. The value of d n can easily be calculated without resorting to the 
first proof: a straightforward inspection shows that the expressions (2.5) and (2.9) 
both induce an asymptotics of the form 

D(<Ji, .. ., <7 n _i, cr n ) ~ n n <j™~ 1 D(<Ji,... ,<j„_i) (a n -> 00); 

the hrst one gives K n = 2, the second one K n = d n /d n -\. From D(a 1) = 2 we thus 
get d n = 2". 


3. Singular Values of Bordered (Skew-)Symmetric Matrices 

In preparation for what follows, in this section we study an algebraic device that 
allows us to untangle the interlacing (2.4b) of even and odd singular values, namely 
bordering a skew-symmetric matrix A £ R” xn with a column vector: A 1 —> (b A). 
By looking at the purely imaginary Hermitian matrix iA we see that each non-zero 
singular value of A occurs with even multiplicity. That is, with n = 2 m + /i, the n 
singular values of A can be arranged as the sequence Si, Si, s 2 , s 2 , • ■ •, s m , s m and, 
if y = 1, also s m +1 = 0, decreasingly ordered according to 
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The results of this section are twofold. First, Lemma 3.1 shows that the singular 
values of {b A) are of the form (with the value = 0 only formally added to the 
list of inequalities if p, = 1) 

tl > Si ^ t-2 > S2 > • • • ^ tm ^ Srh- (2.4b) 

That is, bordering A by a column b splits the double listed pairs ( Sj,Sj ) of singular 
values into ( Sj,tj ) and, if /i = 1, modifies the surplus singular value Sm = 0 into 
some tm subject to the interlacing (2.4b). They are strictly interlacing if 

t\ 5 Si > t‘2 iX £*2 ^ 1 “ tm 5 Srh- 

Second, Lemma 3.2 establishes a coordinate change ( t , s) K > (r, s) through an explicit 
algebraic map such that strict interlacing of t with s corresponds to strict positivity 
of the components of r. 

To begin with, there are orthogonal matrices U and V such that (the last row 
and column of the block partitioning are understood to be /x-dimensional, meaning 
that they are missing if p = 0) 


(S 0 0\ 

UAV' = 0 S 0 L 

\0 0 0/ 


U(b A) 



(uS 0 0 \ 
u 0 S 0 , 
\?7 0 0 0 / 


(3.1) 


with S = diag(si,..., s m ) built from the singular values of A and the partitioning 


( U \ 

ub= U , M,»er, w. 

Hence, the singular values of (b A) are given by the following lemma. 

Lemma 3.1. Let S = diag(si,..., s m ) be a diagonal matrix and u,v £ R m , rj £ 

(fj, = 0,1). Then, with to = m + p, the singular values of the (to + to) x (to + to + 1) 
block matrix 

( a 5 0 0\ 

U 0 50 
\t] 0 0 0/ 

are s i,..., s m , t \,..., tm, satisfying the interlacing property (2.4b). Here, the tj are 
the singular values of the mx (to+1) bordered matrix (r S) with S = diag(si,..., Sm) 
and rj = -fuj + Vj (j = 1,, m); if fi = 1, then r m+ i = |iy| and s m +i = 0. Further, 
there holds 

tl + ■ ■ ■ + t% = r\ H + + si + ■ ■ ■ + s^ (3.2) 

and, if p = 1, 


tl tm rm ‘ Si ' ‘ - S m - 


( 3 . 3 ) 
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Proof. It suffices to prove that the two matrices 

(uS 0 0\ / 0 S 0 0\ 

= J v 0 S 0 J , M 2 = J r 0 S' 0 I 
\r? 0 0 0/ \\r}\ 0 0 0/ 

with r.j = ^Juj +v?, j = 1,..j m, have the same singular values. Using a Givens 
rotation Uj with 



one gets 

rr ( U 3 S 0 0 W 1 0 \ = / 0 Sj 0 \ 

' V.i 0 sj {0 Uj) Vo 0 Sj) ■ 

Hence, by successively applying these two-dimensional orthogonal operations to the 
corresponding rows and columns (and addressing a possible sign change of the last 
row if n = 1) one transforms M\ into M 2 while leaving the singular values invariant. 

We note that the decreasingly ordered singular values Si,... ,s„ of a matrix 
S and those of the bordered matrix (r S), ti,...,tm, are generally known [13, 
Cor. 7.3.6] to be interlacing as in (2.4b). 

To finish, (3.2) follows from expressing the Frobenius norm of (r S) in terms 
of its singular values t and (3.3) follows from doing the same, if /i = 1, for the 
magnitude of the determinant of that matrix with the last column (which is all zeros 
then) deleted. □ 

The next lemma shows that one can uniquely solve the inverse problem t, 1 —> r 
for strict interlacing. 

Lemma 3.2. With the notation as in Lemma 3.1, let si > S 2 > • • • > s^h ^ 0 with 
Sfh = 0 if n — 1. and let S = diag(si,..., s*)- Then, the map 

$ : r 1 —> t = the decreasingly ordered singular values of (r S) 

defines a diffeomorphism 

$ : R* —► {t G R* : t is strictly interlacing with s} . 

If t is strictly interlacing with s, its preimage r = < E>~ 1 (t) is the unique positive 
solution of the system 

rh 2 

Y) +2 [_ fc 2 = 1 (j = l,...,m), (3.4) 

k=l J 

which is explicitly solved by 
_wt(s|) 


w a (0 = (t-st) ■ ■ ■ (£-4), <*(0 = (£-«?) • • • (£-&)• (3.5) 
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The Jacobian of the inverse map $ 1 is given by 

det (= 1 . A (4»-••»*!) 

\ dt k) l^j,k^rh n-'-Tm (Sl"-S m )^ A( S ^,..., S 2) 


(A* = 0,1). (3.6) 


Proof. The squares of the singular values ti, ..., of (r S) are the eigenvalues of 

(r S)(r S)' = SS' + rr' = diag(sf,..., 4) + rr'. 

Now, any set of values t? for which tj satishes the interlacing property (2.4b) can be 
obtained in this way, that is, as the eigenvalues of a positive semi-definite rank-one 
perturbation of D = diag(sf,..., 4) (see, e.g., [18, Sect. 2]). Since rr' does not 
depend on the signs of the individual entries of r, we can always choose r G R™ 0 . 
If r„ = 0 for some v, then the ^-th row and the u-th column of rr' are zero which 
means that sf appears among the values of tj. Hence, strict interlacing implies 
r G M> 0 - 

Given such an r G R™ 0 , th e eigenvalues t? of D+rr' are known ([12, Lemma 8.4.3]) 
to be strictly interlacing with the and satisfy the secular equation 

f{tj)= 0 (j = l,...,m), /(A) = l + r'(D- A/) _ 1 r, 

which is (3.4). Since the determinant (3.7) given below is non-zero and, hence, the 
Cauchy matrix 



is non-singular, there is a one-to-one correspondence of r G with those t that 
strictly interlace with s. Because each of the steps 1 K > C H► r is smooth, we have 
therefore proved that is a diffeomorphism. 

By relating Cauchy matrices with Lagrangian polynomial interpolation, Schechter 
[16, Eq. (16)] gave a short and simple proof of the explicit formula (3.5). Differentia¬ 
tion with respect to tk gives 

t dr 3 _ r jtk 

lk ~ at, ~ «; - 1\ ■ 

Hence, 


J = diag(ri,.. .,rm)C diag(H, • • 


which implies det J = t\ ■ ■ ■ f* rq ■ • • r„ det C. Now, using the explicit determinantal 
formula [16, Eq. (4)] 


det C = 




n 0 vq-si) 


(3.7) 


= (-l) 


m Utjsl) ■ ■ ■ UtjSfo) 

4(4 ••• 4 ( 4 )’ 


and 
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together with the following straightforward evaluations of the product terms 
(-1 = Wit) - s 2 k ), w'(s?) ■ ■ • w'(s£j = ]J(sf - 4 )t 


one gets 


With 




, j _ t\ ■ ■ -tfh Ylj< k ( t2 j -tl) 
" 4) 


~ t i)= 

j<k 




I]>j “ S k) = A ( s2 rh, ■ ■ ■ , 4) = (si---S m ) 2/i A(> 
j<k 

one finally gets the expression (3.6) by using (3.3) if p = 1- 


mvrl 


), 


□ 


4. Random Matrix Models for the Odd and Even Singular Values 

Because of interlacing, the factorization of the joint density stated in Theorem 2.1 
does not reveal an independence between the x and the y components of the singular 
values, or to the same end, between the t and the s components. If we change, 
however, the (t, s) coordinates to the (r, s) coordinates introduced in Lemma 3.2, 
the interlacing is replaced by just a positivity condition on the r components. The 
following theorem, which sharpens Theorem 1.1, shows that not only are the r and 
the s components independent of each other but both sets of components have so 
much additional structure that they can be completely described in terms of known 
distributions. 

Theorem 4.1. Applying the transform (f, s) i —> (r, s) of Lemma 3.2 to the (t, s) 
parametrization (2.4) of the decimated ensembles odd|GOE n | and even|GOE„| 
defines a set of random variables r k , which are distributed as \2 for k = 1 ,..., to 
and, if pL = 1, distributed as for k = m + 1. They are independent of each other 
and of the even singular values s, which are jointly distributed as 

even |GOE„| = aGUE n . 


Proof. In terms of the (t, s) coordinates, the joint density (2.6) of the singular 
values of the GOE can be recast in the form 


q(s;t) 


c n 2 n n\ • (si • • 


SmY {ti ■ ■ -trh) 1 M 
■A(s 2 m ,...,s 2 1 )A(tl,...,t 2 1 ) 


2 .2 
g j — 1 2 j = l 2 


(4.1) 


where the case distinction between even (/z = 0) and odd (/z = 1) orders n has been 
expressed in terms of powers. If we apply the coordinate change ( t , s ) i —> (r, s) of 
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Lemma 3.2, which is a diffeomorphism up to an exceptional set of zero probability, 
the density with respect to the (r, s) is 




• q(s; t) = n • • • r m 


{sx-'-SmY 



m \ 

5^c n 2 n u\ • JJ • A {s 2 m ,..., sj) 2 j 


with Sfj, = (7 t/2 ) m ^ 2 . Here we used expression (3.6) for the Jacobian and simplified 
the exponential functions according to (3.2). On their supporting domains, the first 
m factors of the resulting density are a % 2 -density each, the next one is a Xi-density 
if p = 1 (disappearing if p = 0), and the last one is the joint density of the anti-GUE 
of order n, see [15, Sect. 13.1] or [9, Ex. 1.3.5(iv)]. □ 


Remark 4.1. As a side product, the proof shows that the normalization constant 
a n of the joint density of the anti-GUE, if extended by symmetry to be supported 
on [0,oo) m , is given by 


(t2 n n\ 
a n — c n ^ J ^T 

This is consistent with the explicit 
Eq. (1.163)/Eq. (4.157)]. 


(n = 2 m + n, n = 0,1). 
formulae for c” 1 and a” 1 


given in 


[9, 


The proof of Theorem 4.1 shows that the joint density p{t\s) of the t variables 
conditioned on s is given by the expression 

'‘■"A (t 2 


p{t\s) = 


1 (fr-dm) 

Sfi (Sl • • • SmY A(S^ 


A 

v m ’ ‘ 


di) 

,s'i) e 


= u)/2+Er=i u/2 


(4.2) 


This is just a particular case of a general result by Forrester and Rains [11, Cor. 3], 
which gives the probability p(t\s) if the tk are the solutions of the secular equa¬ 
tion (3.4) with parameters r 3 being independently gamma distributed. In retrospect, 
we could thus have proved Theorem 4.1, based on Theorem 2.1, starting with the 
Forrester-Rains formula (4.2) and working backwards. 

Now, Theorem 4.1 and Lemma 3.1 yield a new random matrix model for |GOE„| 
which amounts to the singular values of certain bordered skew-symmetric Gaussian 
matrices. 


Corollary 4.1. Let X £ R raxn be a random matrix with independent standard 
normal entries. Denote by G = (X + X')/2 its symmetric and by A = (X — X')/2 its 
skew-symmetric part. Let T n be a Xn~distributed random variable that is independent 
of X. Then, both the singular values of G and the singular values of the bordered 
matrix (with e\ denoting the first unit vector) 

H = (r„ei A) 


(4.3) 



14 F. Bornemann and M. La Croix 


are jointly distributed as those of the GOE of order n. The same holds if the matrix 
H is obtained from bordering A with an independent standard normal vector. 

Proof. Note that the singular values of the symmetric part G = (X + X')/2 are 
by definition jointly distributed as the singular values of the GOE of order n. 

As discussed in the derivation of (3.1), the singular value decomposition of A 
takes the form 

(S 0 0\ 

UAV’ = | 0 S 0 J , S = diag(si,. .. ,s m ), 

\0 0 0/ 

where the last row and columns are missing if n is odd. By symmetry, the orthogonal 
matrix U is Haar distributed—independently of {si,..., s m }, which are jointly 
distributed as the anti-GUE of order n, cf. [15, Sect. 13.1] or [9, Ex. 1.3.5(iv)]. 
Applying U to the first column H\ of H defines 

( U \ 

UHi = T n Ui = I v I , u, v € R m , ry € 

Since the first column U\ of U is uniformly distributed on the sphere S'" -1 and r„ 
is independently ^-distributed, we see that UH i is a standard normal vector, see, 
e.g., [5, Sect. V.4]; the same conclusion holds for standard normal H 1 . Hence, the 
variables 

Vj = y Cf+v] = 

are independently X 2 -distributed and, if p = 1, r m + 1 = |?y| is independently %i- 
distributed. Comparing the results of Lemma 3.1 and of Theorem 4.1 finishes the 
proof. □ 


The random matrix model of the last corollary can easily be turned into the 
following sparse model that separates the even and odd singular values. Note that 
just one of the two matrices is square, the other is rectangular. 


Corollary 4.2. Let n = 2m + p, p = 0,1 and let t±, ... ,t„ be independent random 
variables, with t & distributed as Xk ■ The union of the singular values of both the 
bidiagonal matrix 

B° dd = ( Tn e i H* ven ) e R (™+M)x(m+i) 
and the bidiagonal matrix H® ven € R( m +A xm ; defined by 


B 


even 

0 


l 

7^ 


/ 72m— 1 
72m—2 72m-3 




V 


7-2 T\J 


^even 


/ 7~2m 


\ 


l 

7^ 


72m—1 72m—2 

r 3 t 2 

tJ 


\ 
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is jointly distributed as |GOE„|. Here, the singidar values of B° dd correspond to 
odd |GOE„| and the singidar values of J3® ven correspond to even |GOE„|, both drawn 
from the same ensemble. 


Proof. Using the notation of Corollary 4.1, a Householder tridiagonalization of A, 
rescaling rows and columns by —1 as necessary, yields UAV' = T with orthogonal 
matrices U, V and 


/ 0 r„_ i 

I T~n— 1 0 T n _2 


T = 


1 

V5 


Tn —2 0 7"n—3 




j t -2 On 

V -1 0/ 


where the entries are jointly distributed as independent y^-variables with degrees 
of freedom k = 1,... , n, see [6, Sect. II]. Since Householder triadiagonalizations do 
not operate on the first row and column, we have Ue i = e\ and, hence, 


U ■ (r n e i A) ■ 



(r„ei T). 


Thus the matrices (r n e i A) and (r n ei T) have the same singular values. By a 
simultaneous row and column permutation of T so that the odd columns and rows 
occur before the even ones, we see that the matrices [r n e\ T) and (with the length 
of the first unit vector e\ adjusted) 

/r„ ei 0 B- en \ 

V 0 (B^y 0 ) 

have the same singular values. Interlacing shows that the singular values of B® ven 
correspond to the even ones of (r„ei A) and the singular values of (r„ei H® ven ) 
correspond to the odd ones. □ 


Remark 4.2. The mapping from the singular values of a sample drawn from the 
GOE to r coordinates is deterministic and can be made explicit (the same remark 
applies to the construction of the £ variables in the next section): starting with (r, s) 
coordinates, the r are obtained by applying appropriate orthogonal row and column 
transformations to the matrix 

h \ 

Sm 

o ••• o / 


0 0 S' 
r-S 0 


where S = 


if p, = 0, the last row of S is missing. 
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5. Square Bidiagonal Matrix Models and the Determinant 

To study the distribution of determinants we turn the bidiagonal random matrix 
model of Corollary 4.2 into one with square matrices only. Key to this transformation 
is the following variant of a result by Dumitriu and Forrester [6, Claim 6.5]. 


Lemma 5.1. Let the variables Tk (k = 1,..., 2 m — 1) be distributed as Xk, with the 
distribution of 72m arbitrary such that n,..., T 2 m are independent of each other. 
Then the singidar values of the m x (m +1) bidiagonal matrix 


B = 


! 72m T~2m -1 

72m—2 72m—3 


\ 7-2 tJ 

are the same as those of the m x m bidiagonal matrix 

^2m+l f ) 2m—2 ^ 

^2m-l f,2m—4 


R 


& ?2 

u 


(5.1) 


(5.2) 


constructed by the normalized reduced RQ-decompositionf B = RQ with a row- 
orthogonal matrix Q, that is, by the almost surely positive solution of the set of 
equations 

^ 2 k+i + ^ 2 k -2 = T 2 k + T ik-i (k = l,..., m), (5.3a) 

6fc+i6fc = 72fc+iT2fc (k = 1,..., m - 1). (5.3b) 


The variables £ 2 , • • • ,f, 2 m-i are distributed as X 2 , ■ ■ ■ , X 2 m-i/ they are independent 
of each other and of T 2 m - The variable f, 2 m+i is of the form 

6m +1 = + where £1 = tv^ - & m _ 2 

is distributed as Xi and is also independent of f 2 , ■ ■ ■, ^ 2 m-i and of T 2 m - 


Proof. A well-known result [5, Thm. IX.3.1] about the ^-distribution states that 
the involution 

<t>(X,Y,Z)=(z Jc J ^ y ,Z^ y ,X + Y^j (5.4) 

maps a set of mutually independent random variables X, Y , Z distributed as xf, 
Xs and Xr+s 1° a new se t °f mutually independent random variables of exactly the 


b The /^-factor can equivalently be obtained from the Cholesky-type decomposition RR' = BB'. 
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same type. Starting with rip = ti, the system (5.3) is recursively solved for the 
variables 6, • • •, 6m -1 by 

(rlk+i^2k^2k+i) = H T lk, T 2 k ,T% k+1 ), (k = 1,... ,m - 1). 

Hence, the variable 6 = Ti jm and the thus constructed 6, • • • , 6m-i are independent 
of each other and of the not yet used variable T 2m ; they are distributed as \k 
(k = 1,..., 2m — 1). Because <j> is an involution, there is 

— T 2 - 4 - f 2 — f 2 A- £ 2 

Hence, the yet to be used k = m case of Eq. (5.3a) finally implies the asserted form 

of6m+i- □ 


Remark 5.1. The use of the involution (5.4) has been motivated by the observation 
[7, Lemma 1] that for 2x2 matrices the i?.-factor of the RQ-decomposition of a 
lower triangular matrix 


L = 



= RQ, 



is induced by the transformation (£ 2 , ry 2 , f 2 ) = (f{x 2 , y 2 , z 2 ). Basically, the 17-factor 
(5.2) of the bidiagonal matrix (5.1) is then obtained by successively applying this 
transformation along the diagonal from the lower right to the upper left. 


Application of this lemma to Corollary 4.2 yields sparse random matrix models 
for the odd and even singular values of the GOE in terms of bidiagonal square 
matrices. We begin with the case of even order n = 2m. 


Theorem 5.1. Let 6, • • • ,6m be independent random variables, with 6 distributed 
as Xk ■ The union of the singular values of the two bidiagonal square matrices 


and 


oodd 1 

= 71 


( \ A 1 -C... &m-2 \ 

6m—1 6m—4 

6 6 

6/ 


neven 

R ° = 71 


(£l (,2m —2 ^ 

(,2m—1 (,2m —4 


^5 (,2 

6/ 


is jointly distributed as |GOE 2m |. Here, the singular values o/l?g dd correspond to 
odd|GOE 2m | and the singular values of R^ ven correspond to even|GOE 2m |, both 
drawn from the same ensemble. 
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Proof. Let ti,. .. ,72m be independent random variables, with Tk distributed as 
Xk- Prepending a zero column to the second matrix, Corollary 4.2 shows that the 
singular values of the two matrices 


B 


odd 

0 


l 

71 


/_odd _ 

' r 2 m r 2m-l 


72m—2 72m—3 


\ 


V 


72 Tl/ 


and 


where 


/,_even _ 

T 2m T 2m-1 


neven 

" 71 


72m—2 72m—3 


7 2 7i/ 


72m d = 


even _ 
'2m — 


are jointly distributed as |GOE 2m |- Here, the singular values of Hg dd correspond to 
odd|GOE 2 m| and the singular values of B g ven correspond to even|GOE 2 m|, both 
drawn from the same ensemble. If we apply the construction of Lemma 5.1 to both 
matrices Hg dd and Rg ven simultaneously, we obtain the f?.-factors f?g dd and i?g ven 
with one and the same set of variables £i,..., £2 m subject to the asserted properties 
and additionally 

£ 2 °m + l = V '(! + (7 2 ° dd ) 2 = ^ e2m+l = ^ = £l • 

By defining f 2 m = 72m we thus get the asserted form of f?g dd and i?g ven . □ 


Remark 5.2. The matrices i?3 ven of Theorem 5.1 and Rg ven of Corollary 4.2 are 
superficially related, up to a different sample of the independent random variables, 
by a transposition followed by a cyclic permutation of their diagonals. Such a 
transformation would not, in general, preserve singular values, but depends instead 
on properties of the x distributions. 

The case of odd order exhibits a similar structure. The equivalence between 
the models Bf ven and i?| ven used in the following theorem is also noted, for the 
anti-GUE, in [7, Sect. 2] and [6, Clm. 6.5]. 


Theorem 5.2. Let £i,...,^ 2 m+i be independent random variables, with £*, dis¬ 
tributed as Xk ■ The union of the singular values of the two bidiagonal square matrices 


j^odd 


l 

71 


/ 72£l 72^2 m 

6m+l £2 m —2 


\ 


V 


£5 £21 
£3/ 
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and 


^>even 


1 

71 


/6m+l 6 m —2 

^2m-l ^2m-4 


V 


£5 


\ 


6 I 

6/ 


is jointly distributed as |GOE 2m+1 |. Here, the singular values of R° dd correspond 
to odd|GOE 2m+1 | and the singular values of R\ ven correspond to even |GOE 2m+1 |, 
both drawn from the same ensemble. 


Proof. Let t 15 . .. ,r 2m+1 be independent random variables, with t*, distributed as 
Xk- By transposing both matrices and prepending a zero column to the first one, 
Corollary 4.2 shows that the singular values of the two matrices 


B 


odd 

1 


1 

71 


/0 v / 2r 2m +1 

i t 2 m T 2m-1 


\ 


V 


T 2 Tl) 


and 


^even _ 


72 


! T 2 m 1~2m— 1 

T~2m-2 T2m-3 

V T 2 Tl) 


are jointly distributed as |GOE 2rn+ i|. Here, the singular values of Hg dd correspond 
to odd |GOE 2m+ i| and the singular values of Bg ven correspond to even |GOE 2m+ i|, 
both drawn from the same ensemble. If we apply the construction of Lemma 5.1 
simultaneously to Bf ven and to 


^odd 


1 

71 


/_odd _ 
r 2m+2 r 2m+l 

T2r 


T2m-1 


\ 


V 


T2 Tl) 


odd 

'2m+2 


= 0 , 


which is B° dd with its first row rescaled, we obtain the .R-factors I?® ven and 


jj^odd 


l 

71 


/&m+3 6m \ 

6m+1 6 m —2 1 


l 6 )) 
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with one and the same set of variables £ 1 ,..., £ 2 m +i subject to the asserted properties 
and additionally 

^ d + 3 = ^? + « d + 2) 2 =^l- 

By restoring the proper scaling of the first row, we thus get the asserted form of 
R\ ven and Rf d . □ 

Remark 5.3. This theorem implies that, for odd order n = 2m + 1, the product 
det R° dd of the odd singular values of GOE„ and the product det i?f ven of the even 
ones are related by 

det Rf d = Ci det i?r n , 

where Ci is a random variable distributed as Xi which is independent of even |GOE n |. 
This is nothing but (3.3), recalling that by the proof of Theorem 4.1 the variable 
r m+ i is distributed as %i- 

As an immediate consequence of the preceding theorems, stated in the following 
corollary, the magnitude of the determinant of the GOE can be expressed as a 
product of independent random variables. Even though Tao and Vu speculated that 
such a representation does not seem to be possible [17, p. 78], a precursor of this 
result was recognized implicitly by Delannay and Le Caer who noted in [4, p. 1531] 
that the Meijer G-function they used to describe the distribution of the determinant 
when n is odd could be sampled as a product of independent gamma-distributed 
random variables. They did not, however, have any interpretation for these variables 
in terms of the underlying ensemble, nor did they recognize the possibility of sampling 
when n is even. 

Corollary 5.1. Let G n be drawn from the GOE of order n = 2m + p with parity 
p = 0,1 ,rh = m + p. Then the determinant of M n = \/2G n factors into independent 
random variables of the form 

\detM n \ = • £3 • £5 • • • C 2 A-I! (5-5) 

with 

1 i 1) = £1 • + 2£2 [p = 0), J?i 1) = V^£i ( 11 = 1). 

Here, the £*, are mutually independent random variables distributed as Xk- 

Proof. The assertion follows from the observation that \detM n \ is the product of 
the singular values of M n and therefore, by Theorem 5.1 and 5.2, distributed as 

det(\/ 2 i?® ven ) • det(v / 2 i?° dd ). 

Multiplication of the diagonal terms of the bidiagonal factors finishes the proof. □ 

In retrospect, once we know that |detGOE„| is distributed as a product of m 
independent random variables, all the factors can be readily identified in the Mellin 
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transforms computed by Delannay and Le Caer in [4, Eqs. (26/41))], although, for 
even n = 2 to, the Mellin transform of the factor 


V2m ~ Cl \j + 2^2m 

into the expression 


n3(s—1)/2 ^ ( 2 ) ^ ( 5 + m 2 ) 

r(i)r(| + m) 


F l 3 1_S - 

2^1 01 o ’ 


S 1 

2 + ’" ; 2 


(5.6) 


with the hypergeometric function 2 F 1 may not be familiar to most observers. This 
aspect of their paper has been missed by several commentators, with Mehta, in [15, 
§26.6], omitting their expression for even n = 2m since the inverse Mellin transform 
of (5.6) cannot be readily written down. 


Remark 5.4. In the case of odd order n = 2 m + 1, the density of det GOE 2m +i is 
necessarily odd, since the eigenvalue density (2.1) is even, but det(G) = — det(— G). 
It follows that in this case the sign of the determinant is statistically independent of 
its magnitude, and we can obtain the distribution of the determinant by replacing 
£1 by a standard normal variable. No corresponding result is available for even order 
n = 2m, although the factored presentation of the odd moments of det GOE 2m in 
[1, Eq. (23)] suggests that the distribution of the determinant should involve many 
of the same factors. 


6. Central Limit Theorem for the Determinant 

Delannay and Le Caer used an explicit computation of the Mellin transform of the 
even part of the distribution of det GOE n to derive the cumulants of the potential 
V = log |detGOE„|, and to show that V is asymptotically Gaussian [4, Section III]. 
Tao and Vu extended this log-normality to determinants of a wider class of Wigner 
matrices, and provided an alternate proof in the case of Gaussian matrices in [17]: 
based on analyzing tridiagonal sparse models for the GOE and GUE eigenvalues, 
they found a way to approximate the log-determinant as a sum of weakly dependent 
terms, which then yields the asymptotic log-normality by stochastic calculus and 
the martingale central limit theorem. In this section we present yet another, much 
simpler proof based on the factorization of the magnitude of the determinant into 
independent random variables. In particular, our proof elucidates the difference 
between the GOE and the GUE in the scaling of the central limit theorem. 

We start by recalling that, parallel to the factorization given in Corollary 5.1 
for the GOE (/? = 1), Edelman and La Croix [7, Tlrm. 2] obtained a factorization 
for the GUE (/3 = 2): with G n drawn from the GUE of order n, the determinant of 
M n = \[2G n factors into independent random variables of the form 

|det M n | = 7 • £ 3 £ 3 • £ 5^5 • • • ^ 2 m— 1 ^ 2 m —1 (6.1) 

with 

Vn ] = £i£n+! (m = 0), ?^ 2) = Cl (/* = !)• 
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Here £ 1 ,... £ 3 ,..., £, 2 m-i are mutually independent random variables with both 

and fk being distributed as Xfc ■ Note that, except for the (asymptotically irrelevant) 
change in the factor rj^\ the transition from the GOE to the GUE just amounts for 
splitting the terms into the products of independent factors. It is precisely 
this split which causes the appearance of [3 in the denominator of the central limit 
theorem when written in the following form. 


Theorem 6.1 (Tao and Vu [17, Thm. 4]). 

holds, as n 00 , the central limit theorem 

log | det M n | - § log re! + ^ log n d 


With the notation as above there 

N( 0,1) (/3 = 1,2), (6.2) 


where —> denotes convergence in distribution. 


To prove this theorem for f3 = 1 and (3 = 2 in parallel, we split 
log|detM n |=y n <«+Z(« (/? = 1 , 2 ) 
into the random variables Yn^ = log r 7 ^ and Z^ dehned by 
Z n^ = 2 log ^3 + 2 log £5 + ••• + 2 log f, 2 rh— 1 > 


Z n ] = ( lo § 6 + log I 3 ) + (log £5 + log | 5 ) + • • • + (log &rh -1 + log l) , 

We immediately observe the relations 

E (z* 1 *) = E (z< 2 >) , Var (z«) = 2 Var [z^) . (6.3) 

Note that the factor of two between the variances is caused, in the transition from 
GOE to GUE, by the above mentioned split of £ 2 into the product fkfik- 

Now, while proving the central limit theorem in the /? = 2 case, Edelman and 
La Croix [7, Cor. 2] obtained, in passing, the following result. 


Lemma 6.1. The random variable Z^P satisfies, as n —► oo, a central limit theorem 
of the form 



Zn ] ~ \ log n\ + V log n+\ log n ^ 

y/j> lo & n 


(/? = 1,2). (6.4) 


Proof. The proof of [7, Cor. 2] proceeds, first, by establishing asymptotic expansions 
based on explicit calculations of the mean and variance of log x-distributed variables, 
namely, 

E ( z i 2) ) = \ log( 2 ?h - 1 )! - Cjf l og ( 2 m - 1 ) - * log( 2 ?h - 1 ) + 0 ( 1 ), 
Var (-Zd 2) ) = \ log(2?h - 1) + 0(1), 
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( 2 ) 

and, next, by showing that Z„ ; satisfies a Lyapunov condition of order four. Hence, 
the Lindeberg-Feller central limit theorem can then be applied to Z\ ' and gives, 
by noting that 

log(2m — 1)! = log ?x! — (1 —/z)logn, log(2m — 1) = logn + 0(1), 

the asserted limit (6.4). By realizing that the sums and Z^ basically share 
the same Lyapunov condition, the central limit theorem for zV can be induced 
from that of Z\ ; by means of (6.3). □ 

The difference between the central limit theorems of log det M n | and of zV 
enjoys the following strong convergence result. 

Lemma 6.2. The random variable yV' 1 satisfies, as n —> oo, 

y(£) __ hzM I can 

= " r-^ g (M = 0,1), (6.5) 

V/3 lo s n 

where denotes almost sure convergence. 

Proof. The case /.t = 1 is trivial, since in that case yV is independent of n. Applied 
to a sum of squares of independent standard Gaussians, the strong law of large 
numbers gives, as n —> oo, 

nr i£ E ^ = 1; andi hence, n _1 (£\ + 2^) 2. 

Taking the logarithm gives 

log£n - ^logn -^4 0, log^7^-ilogn^ilog2, 
which implies the assertion for /x = 0. □ 


Now, adding (6.4) and (6.5) gives, by Slutsky’s theorem, 
log |det M„| — pog n! + ^ log n = + ^ ^ 

Y lo g n 

which finishes the proof of the central limit theorem (6.2). 


(0 = 1 , 2 ), 


7. Integrating Out the Odd or Even Singular Values 

Here, we present another proof of Theorem 1.1. If we are interested only in the 
distribution of the even singular values, then it is possible to proceed from the joint 
probability density (4.1) by integrating out the odd ones. While not exposing any 
additional structure, such as in Theorem 4.1, this approach is conceptually more 
straight forward, and offers the advantage that it can also be used to establish the 
determinantal formula (7.4) for the probability density of the odd singular values. 
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This is of interest in its own right, since we constructed separate sparse random 
matrix models for the odd singular values in Corollary 4.2 and in Theorems 5.1 and 
5.2. Moreover, the technique extends to the symmetric Jacobi and to the Cauchy 
ensembles [3]. 


7.1. Integrating out the odd singular values 

Recalling (1.5), we rewrite the expression (4.1) of the joint density in the form 
q(s; t) = c n 2 n n\ • ^ M (si,..., s m ) ■ gi-^h, ...,**) 

with functions 

m 

g a (zi,-..,z m ) = z%e~ z * 12 ■ (7.1) 

fc=l 

Corollary 7.1 below shows that integrating out the odd singular values t subject 
to the interlacing (2.4b) gives the following marginal density of the even singular 
values with 5 M = (7r/2)^/ 2 : 

^even(si, * • • , Sm) — Sfj,C n 2 77.! • ^(.S1, . . . , S m ) 


= 5^c n 2 n n\ ■ s^e s * ■ A(s^,..., s?) 2 . (7.2) 

k =1 

Since the last expression is the joint density of the anti-GUE of order n, see [15, 
Sect. 13.1] or [9, Ex. 1.3.5(iv)], this is nothing but Theorem 1.1 spelled out in terms 
of densities. 

The integration is based on the following lemma and its first Corollary 7.1. 


Lemma 7.1. Let 


( x^e x2 / 2 \ 


£°(*) = 


x K+2 e -x 2 /2 


r K+2n—2 p —x 2 /2 




(«= - 1 , 0 , 1 ), 


with the understanding that, instead of x x e x2 / 2 7 the first entry of e^(x) is the 
expression 

\V2 

Then, for n = 0,1, there holds the integration formula 

C X n +1 


V-i( x ) = -d 2 erf 


r X 2 rXn + l , v 

/ d£ n det(4 n) (£i) ... e W(£ n )) 

j Xl J Xn ' ' 


••• ei-lf^n+lA 
, 1 ••• 1 ) 


= det 
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Proof. Integration by parts yields the three-term recurrence of antiderivatives 

J e~^ /2 df = - 77 -iOr), 

J X e +1 e~ e/2 df = + k j" e-'e-S 2 ' 2 d£ (, k = 0 , 1 , 2 ,...), 

and, hence, by simplifying notation to e K (x) = eV(x), 

J e K (£)ci£ = L K e K _iW (k = 0,1) 

with a lower triangidar matrix L K G R nx ™ having —1 all along its main diagonal. 
We thus calculate 

/ %2 PXn+1 

d £!••• / df n det(e K (^) ••• e K (£ n )) 


= det / e K (£i)d£i •• • 




/ X 2 

Xn -\-1 \ 

e K -i 

•• e K _i 

\ Xl 



= (-!)" 


X2 

Xn-\-1 \ 

e K-l 

■ e K _i \ 

Xl 


0 

0 / 


e K _i(xi) e K - 1 (x 2 ) ••• e K _i(a: n+ i) 


In the last step we added the first column to the second, then the second to the 
third, etc. □ 


Corollary 7.1. Let g^ be as in (7.1) and put Sm = 0 if n = 1. Then, one has the 
integration formula 


fOO f-S 1 


/ dt i dt 2 •••/ dtrhg 1 ^^{t 1 ,...,trh)=5^g^{s 1 ,...,Sm) (m = 0 , 1 ) 

J Si J S 2 J Sm 

with 6 n = ( 7 t / 2 Y^ 2 . 

Proof. Using the notation of Lennna 7.1, we first observe that 

9 ^z u ...,z m ) = det (e^ (z m ) ••• e^ m) (. 2 i)) . (7.3) 
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Now, Lemma 7.1 yields, first using e^^oo) = 0, that for g = 0 



and then, using e^ +1 ' l (0) = 0 and e^ +1 ^(oo) = — {y/irj2, 0)', that for /r = 1 


/»00 nS i pSm 

J dti J dt 2 ---J dt m+1 det [e i 0 rn+1) (t m+1 ) ••• e[> m+1) (ti)J 
= det ^ e<T 1 ) («m) • • • ^(-i) e^ + 1 ) (oo) j 


(—!) m det (e5 +1) (a m ) ••• e^ +1) (*i) e^ +1) (°°)) 


( 77 _ i ( s to ) 
e[ m \s m ) 


V-i(si) — x/¥ 
ei m) ( S i) 0 


^det (eC(sm) ••• e[ m) (si)) 


which finishes the proof of the corollary. 


□ 


7.2. Integrating out the even singular values 

The following second corollary of Lemma 7.1 will allow us to integrate out the even 
singular values from the density q(s; t). 


Corollary 7.2. Let g^ be as in (7.1) and put t m +i =0 if g = 0. Then, for g = 0,1, 
one has the integration formula 



nt2 


ds 2 


>t 3 


>tm +1 


ds m g^(si ,.. 


s m ) = det 


MW) 

V ^1—/x(^ra) 


5i-n{ti) J 


with Soft) = 1 and Si(t) = \J 7r/2erf(f/-\/2). 


Proof. Using (7.3) and Lemma 7.1 we obtain 
I ds\ • • • j dsm g^(si ,..., s m ) 

a t 2 ” tm -(-1 


fii 

f tm / 

/ ds i•• 

■ / ds rn det ( 

)t 2 

" tm-\- 1 


• • e 


(m) 


r e^\(t m+1 ) 




= det 
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which is already the assertion for fx — 1. For fi = 0, the assertion follows from further 
calculating 


det 


4 m) (o 


• e«(ti) N 


= det 


'o 4”>( (m ). 



= (-irdet(e«(t m ) ■■■ 4l?(i,)) = (-!)”>det ( ) 

(4 n)\ 

( «1((„) ••• J 

which finishes the proof. □ 


Now, by means of this corollary, the marginal density of the odd singular values 
supported on t\ ^ <2 ^ ^ t-m ^ 0 is given as 


Qodd(tl> • • • 5 tm) — C n 7l\2 • </l—p(tm, . . . , £l) • det 


7 e ( ™ _1) (7 - ) 
c l-u Vm) 


4tyto' 


$1—* * * ^1— 


„ \ (m-i) / 

= C„!2" • det P-' (< ” ) (f ‘ 

7l7l — 


det 


'4”:‘’fe) 


(m-i), y 
e l-u GO 


* * * ^1—p(tl) 


(7.4) 


with 


7mW = t^ + 2 rh_ 2 e _ * 2 / 2 , <y„(i) = 


if /i = 0, 

A) if <‘ = L 


Note that the two determinantal factors differ just in their last rows. It is this 
difference that prevents the expression from becoming a perfect square, which is in 
marked contrast with the marginal density of the even singular values. 


8. Gap Probabilities 

Theorem 1.1 has an interesting implication in terms of gap probabilities, that is, in 
terms of the probabilities 

that the interval J contains exactly k eigenvalues drawn from the random matrix 
ensemble RMT of finite order n, or in some scaling limit. Here, RMT will be the 
GOE, the aGUE or the LUE with parameter a. 

To begin with, by a simple change of coordinates, see [7, p. 8], there holds 

(0, «)) = ERjvfc (0, s 2 )) | a=M _ i (m = 0,1). (8.1) 
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By looking at pairs of consecutive values it is easy to see that the event that exactly 
k values of the decimated ensemble even|GOE n |, n = 2m + p, are contained in 
(0, s ) is given by the union of the events that exactly 2k + p — 1 or that exactly 
2k + p values of |GOE„| are in that interval. Since these two events are mutually 
exclusive and since the singular values of GOE contained in (0, s ) correspond to the 
eigenvalues in (— s, s), we thus get from (1.2) and (8.1) proof of 

4oe*(2* + 0-1; (-*.«)) + 4oe m (2* + 0; ( s , *)) = (o, a )) 

= ^ UE (fc;( 0,s 2 ))| o= ^_i (p = 0,1). (8.2) 

For even order (p = 0), a first proof of this formula was given by Forrester [8, 
Eq. (1.14)]. For odd order (p = 1), Forrester communicated to us further proof of 
the k = 0 case, a remarkable tour de force extending the techniques from [8] based 
on generating functions, Pfaffian calculus, and Fredholm determinants—later he 
was able to use the same approach to establish the general k case; for this and the 
extension to the symmetric Jacobi and to the Cauchy ensembles see [3]. 

We first identified the p = 1 form of (8.2) via a heuristic duality principle based 
on three observations. First, the LUE of order m and parameter a = p — m € N is 
modeled by the eigenvalues of m x m-Wishart matrices W = X'X, where the random 
p x TO-matrices X have independent complex standard normal entries. Second, the 
eigenvalues of W = XX' are those of W padded with a = p — m zeros; that is, the 
(k + o)-th eigenvalue of W is distributed as the k- th eigenvalue of W. Last, since W 
is constructed the same way as W, but with dimension m = m + a and parameter 
a = —a, we are thus led, at least formally, to the duality principle 

+ a; (0,t))| o= _ Q = KWfc; M)la=a • (8-3) 

Extrapolated to general a > —1, it can be taken as a natural definition of an 
otherwise undefined expression. Now, formally evaluating the n = 0 form of (8.2) 
at half-integer values of m and k, and invoking the heuristic duality principle (8.3), 
led us to predict the p = 1 form. Since it held up under numerical scrutiny, trying 
to prove this prediction was a key motivation to our present work. 

As already noted by Forrester [8, Eq. (1.16)], the bulk scaling of GOE and the 
hard-edge scaling of LUE allow us to turn (8.2), as n —► oo, into the limit relation 

-®goe( 2& “1 + 16 (—+ s )) + -^goe( 2fc + At; (— s, s)) 

= E^i{k;(0,TT 2 s 2 ),^- \) (/x = 0,1); (8.4) 

a remarkable formula previously established by Mehta [15, Eqs. (7.5.27/29), 
(20.1.20/21)] using two different, but much more involved methods. In contrast 
to (8.4), which offers many advantages for the numerical calculation of gap prob¬ 
abilities of the GOE in the bulk scaling limit [2, Sect. 5], the finite dimensional 
version (8.2) is not yet a closed recursion that would allow us to calculate the 



The Singular Values of the GOE 29 



Fig. 1. Fluctuation statistics (densities) of 100 000 samples of GOE n , n = 50, vs. the theoretically 
predicted soft-edge scaling limits; top row: largest eigenvalue vs. Fi(l;s) (left), second-largest 
eigenvalue vs. i^fljs) (right); bottom row: largest singular value vs. Fi(l;s) 2 (left); second-largest 
singular value vs. i* 2 (l;s) (right). We have consistently used the 0(n —2 ' 3 )-scaling of Johnstone 
and Ma [14, Thm. 2]. 


gap probabilities of the GOE on symmetric intervals: a complimentary expression 
evaluating 

^GOE^fc - W (-*,*)) + 4oE M (2fc + 1 - W (s, a)) ( M = 0,1) 

is still missing. By the same arguments that justify (8.2) such an expression would 
establish the gap probabilities of the decimation ensemble odd | GOE„ |, whose joint 
distribution is given by (7.4). 

To finish the paper, it is amusing to note that all three cases of the Tracy-Widom 
distributions 

F f5 (l;s)=Ef t (0;(s,oo)) 09 = 1,2,4), 

can be sampled from the soft-edge scaling limit of the spectrum of just the GOE, 
i.e., the /? = 1 case (the case /3 = 2 corresponds to the GUE, /3 = 4 to the GSE), 
see Fig. 1. First, let Ai, A 2 denote the largest and second-largest soft-edge scaled 
eigenvalues of the GOE. In the large-matrix limit, asn-> oo, they are asymptotically 
distributed as 

Ai ~ Ei(l; s), A 2 £f 4 (1;s). 
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The first assertion is the definition of the distribution Fi(l;s), while the second 
follows from the decimation relation, see [10, Thm. 5.2] and [11, p. 44], 

GSE m = evenGOE 2m +i. 

Second, let E l5 E 2 denote the largest and second-largest scaled singular values of 
the GOE. They are asymptotically distributed as 

SE 2 ~ F 2 (l; s). 

Here, the first assertion follows from the asymptotic independence of the extreme 
eigenvalues of the GOE and the second follows from (1.2) as follows: E 2 behaves 
like the largest scaled eigenvalue of the anti-GUE which, like that of the GUE, is 
governed by the Tracy-Widom distribution F 2 (l; s). 
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